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1. Introduction 



The search for the imprints of the string theory in the cosmology is one of the most 
important problems in string cosmology. Cosmic string[|l|, as a promising candidate, has 
drawn much attention in the past few years. The renewed interest |^ in the cosmic string 
resides in the fact that in the low-energy models the cosmic strings could be created after 
inflation. This is in contrast with the prediction of the old Planck-scale models, where the 
cosmic strings produced before the inflation are "blown away" [0] . It has been shown that 
some strings, including fundamental string and solitonic strings, are sufficiently stable and 
a few of them may survive and be observable in near future through gravitational wave 
detectors such as LIGO and LISA. 

In the string-inspired models on inflation, there is a class of models based on the 
brane-antibrane annihilation. In such configurations, the tachyon field develop. As tachyon 
condense, very massive closed string states are created, which could be an origin of the 
cosmic string. However these massive fundamental string states are not stable and may 
decay to other massive or massless states. In order to understand the fate of a typical 
massive F-string state better, it is essential to study its decay systematically. On the other 
hand, even without considering its phenomenological implication, the study of the decay of 
a typical massive string state is an interesting issue on its own right in string theory. From 
the investigation of the rolling tachyon and the geometric tachyon realization in the (NS, 
D)-system0, it has been known that the tachyon condensation will radiate various kinds 
of closed string states [0 , even though there does not exist a well- understood treatment 
of the backreaction yet. To understand the tachyon condensation and the mysterious 
open/close string duality 0, the study of the fate of these massive string states is necessary. 
This is the question we will address in this note. 

The decay of massive strings is an old problem of string theory and it is generally 
difficult to extract detailed knowledge from conventional calculations mainly due to the 
exponentially growing state density. States at the same mass level may be very different, 
especially their decay properties. For example, there exist some peculiar string states, often 
with large angular momentum, which are extraordinarily stable, with lifetime proportional 
to some positive power of their masses [13| [|l^ . 



Fortunately in a realistic situation, such as the rolling tachyon case, details of string 
states are not important, and we are only interested in some sort of averaged decay rates. 
A peculiar long-lived string state may be created but with much smaller probability, it 

1 



must be in an extremal end of a distribution in which the typical string states are amply 
produced and more interesting. The lifetimes of the massive strings are cleverly computed 
in this spirit [|T6[. In [|12[ and |T^, massless string emission from heavy strings were 
considered for both bosonic and superstrings, where the authors average over the initial 
states of the same mass level and fixed four momentum and sum over all possible states 
for the remnant string, a thermal spectrum was obtained. Manes considered in fl^ the 
emission spectrum of the strings of all mass levels for bosonic (open and closed) strings. 
It was found there that the decay rates are universal, in the sense that the leading terms 
depend only on the mass level of the strings involved, and the more detailed properties of 
the states affect only the sub-leading terms. 

To study the stringy properties relevant for cosmology, we need to consider the sit- 
uation where some of the space-time dimensions are compactified. And in the realistic 



models building for early day inflation [|TT| and present day accelerated expansion |jT^ , the 
superstrings rather than the bosonic strings are used. So in this note, we consider first 
the decay process of the bosonic strings with toroidal compactification in Section 2. Here 
the KK modes and winding modes enter and make the process much more complicated. 
And in Section 3, we go on to discuss the superstring process in the flat space where the 



universal expressions similar to [|14[ are found. We do compactification in Section 4 to 
calculate the emission rate in a more realistic string inflation model. Section 5 includes 
our conclusions and discussions. 



2. Bosonic string decay: with compactification 

In this section, we study the emission spectrum of massive closed bosonic strings in 
the case where out of the total D dimensions are compactified on a flat torus. A closed 
string state is characterized by its KK momentum numbers rti, winding numbers Wi, and 
its mass M and momentum k in the noncompact dimensions. The mass shell condition 
reads 

a'M^ = A{Nr - 1) + a'Q\ = 4{Nl - 1) + a'Ql (2.1) 
where Q± = X] ^ "^^^ ? i = 1, ...,dc- For later convenience we also define 

4M|^ = M2-Q|. (2.2) 

The decay amplitude for closed strings can be factorized into the left and the right 
parts, and each part can be calculated in a way similar to open strings. Consider first the 
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string of given mass M, winding Wi and KK momentum rii in its center of mass frame 
decays into two strings with mass Mi and M2, winding wu and W2i, KK momentum riu 
and 77-21, with similar mass shell conditions 



a'M^ = A{NiR - 1) + a'Ql+ = 4(iVii - 1) + a'Ql_, (2.3) 

and 

a' Mi = 4(iV2fl - 1) + a'Ql+ = 4(iV2L - 1) + a'Ql_. (2.4) 
In the following, we set a' — h. Energy conservation gives 



M = ^M^ + k^ + ^jMi + k^, (2.5) 

with k the momentum in the noncompact dimension. 

What we want to consider is the averaged semi-inclusive two-body decay rate. That 
is, for the initial string, we average over all states of some given mass, winding and KK 
momentum. For one of the two final strings, we sum over all states with some given mass, 
winding and KK momentum; only the other string's state is fully specified (by keeping 
explicit its vertex operator). 

This decay rate can be written as 

- ^D-dc 2 Fl Fr j D-3-dc TT D-l C2«^ 
i^serm-zncl- ^2 g{NL) Q{Nr) ^ ^ ' 

with closed string coupling Qc-, compactification radius Ri, and numerical coefficient Ap = 

3-p 



and Fl and Fr are given by 



Fl= Yl Yl \{^f\VL{nu,wu,k)\^i)\^ (2.7) 

Fr= J2 E \{^f\VR{nii,wu,k)\^i)\\ (2.8) 

The vertex operator for the one of the final state strings has been factorized into the left 
part Vl and the right part Vr, whose exponential parts look like e*'^^'^^ , e*'^^'^^, with 
—kj^ = Ml — Qi_, which we define to be 4m|^ for later comparison with open strings. In 

the fiat spacetime the normalization of the initial state and the final state two-body phase 

1,0—3 

space contribute a factor dQ, while the integration over the solid angle Q gives Ad- 



Now momenta in the dc compact dimensions are discrete, so we convert the integration to 
a summation and this results in the factor J^^" R~^. 

The sum in a fixed level is hard to perform, this can be done with a trick: we insert a 
projector Pn over states of level A^, and then sum over states of all levels, in the end we 



obtain a trace |T^] . Let 

dz 



N=f -^^-^ (2.9) 
Ic z 

PnI'^n') = Snn'I'^n'), (2.10) 
so the amplitude squared can be written as 

Fl= (f — z-^^ / ^z'-''^^Tv[z'^vl{k,l)z'^VL{k,l)], (2.11) 
Jc ^ Jc ^ 

where we place the vertex operator at z = 1, takeing advantage of the SL{2, R) invariance 
of the open string tree amplitude. In the following, we define w = zz' , and v = z' . 

2.1. The recursion relation 



Manes observed that the trace in ( p.ll| ) describing a tree level process can be 



carried out by converting it to the two point correlation function on a cylinder with modulus 
t, where 

it = --^. (2.12) 
mw 



Note that 

j d^pTT[z^''V\k,l)z'^''V{k,l)] = f{wf-^{^)^'^ < V^{k,l) V{k,v) >, (2.13) 
where 

oo 

/H = 11(1 (2.14) 

n=l 

The second factor comes from the momentum integral over the zero modes, which does 
not appear in the definition ( |2.7|) and ( p78| ) of Flj^. Remember that compactification does 
not affect the oscillation modes, so we have 

Tr[z^Vl{kLA)z'^V{kLA)] = fiwf-'' < V'[ikLA)Vi^{kL,v) >, (2.15) 
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where the prime means that dJ^Xi^{v) is replaced by v^d!^X>^{l)v~^ in the vertex operator 
so that zero modes are excluded, and D = 26 is the full space-time dimension. Thus the 
master formula for emission of general states including KK modes and winding is 

Fx, = / i d^yN,-N,, < ^/^(^^^ l)V'{kL,v) > . (2.16) 

J W J V 

It can be rewritten as 

Fl=6 —w-'^-awf-^Ir^.-NM (2.17) 





dw 




— w 




w 




dv „ 




V 




V 



V'i{kL,l)V[{kL,v)> (2.18) 

where the contours satisfy \w\ < \v\ < 1, for \z\, \z'\ < 1 . A general vertex operator V{k, z) 
is got from a physical state 

by replacing a'^^ by , and |0; A; > by e*'^^^^), and the normal ordering is assumed. 

We take the gauge* 

^'*^C^l/U2...Mm ~ ^'^^C^lMQ-.-Mm ~ ••• ~ ^'^"^ ^^1 ■ ■ -Mm ~ 0' (2.20) 

for which the contractions of exponentials with dX like terms vanish. Then the correlator 
in Xn{w) factorizes into two parts, the contraction of two exponentials and the contractions 
without exponentials. 

Remember that X scalar correlator on the cylinder is not changed by compactification, 

so 



< a;X^(l)a;"X^(^;) >= -(-)"r7'^^(^;a^,)^+^ In Vi(^;, w), (2.22) 

with 

^ fl - ' 



i;{v, w)^{l-v)ll ^\ ' K (2.23) 



n=l 



* Strictly speaking, such a choice of gauge only exist in the critical dimension: for bosonic 
string D = 26 and for superstring D = 10, which are implied in this paper. 
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In the above correlators, substitution has been made to account for the prime in the vertex 
operator. And dr is vdy for Euchdean proper time t = Inz. The correlator in can thus 
be written as 

< V'[{kL,l)V[{kL,v) >= (2.24) 

with 

n{v,w) =< drXL{l)drXL{v) >= y^n{v^ + ^ ^ ) (2.25) 

1 — 

n=l 



The contour integration in X„ was computed in [jl^. Since the similar calculation is 



to be performed in the superstring case, we in the following repeat the steps of [jl^ . Note 
that 

tjj{w~^v,w) = tlj{v,w), (2.26) 

w 

d:^n{w-\, w) = d!^n{v, w). (2.27) 

As m| = 2{Nl — 1) is even, this makes the minus sign in ( |2.26D irrelevant. The integration 
on the new contour C'^ carries out to be 

i ^y'n < V'[{kL, l)V[{kL, V') >= i / < V'iikL, l)V[{kL, V) > 

which is just 



(2.28) 



J^+^|(«;). (2.29) 
The difference between the two contours is the integration over a small contour around the 



singularity f = 1 or r = Inf = 0. The function '(l){y^ w) has the behavior [|l4| 

tl^iv.w) = -ril + ^r) +0{t^) (2.30) 
near v = 1. So near v = 1, the correlators have the following asymptotic behavior 

< g-*fcLX^(i)g-ifc^x^(«) >=^(y^^)rni = + 1^2^ + o(r2)), (2.31) 

< d-X,il)d-X,iv) >= -dr-'n{vM = + 0(^;^). (2.32) 

Note that one more partial derivative for each vertex operator gives an additional in 
( p.32|) , but at the same time, the corresponding state is more excited, contributing to 
( p.31|) . The two factors exactly cancel 

< V'Hk. l)V'{k, v) >= 4(1 + + 0{l), (2.33) 
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where the coefficient can be fixed to c = 1 by state normalization. So the contour inte- 
gration around f = 1 in our case with compactification is the same as in the noncompact 
case 

/ < v'iikL, l)V[{kL, ^) >^ n + ^. (2.34) 
Jc V 2 

We get the following recursion relation 

2 

2 777 

J,+^2 {w) = w-+^^ [Ir^iw) +n + ^l (2.35) 



which, used in sequence, can be recast into a series summation 

A 

\n — rr?,? (n — ^ 



X^(w) = J][n-mi(p- -)]z/;"P-^"^'(p'-P) + J,(tx;)tx;"^+^"^' (^'+^), (2.36) 
p=i 

where u = n — m\A, with A an integer. Then the contour integration in (|2.16|) can be 
written as 

A 1 1 

F = J](n -ml{p- -))G[Nl - np + -ml{p^ - p)] + F^u. (2.37) 
p=i 

where n = — -/V2L5 and 

FNu=f ^«;-^/(«;)2-^X.(u;)«;"^+^"^'(^'+^), (2.38) 

and the generating function for the mass level density 

00 

Trw^ = f{wf-^ = J2 S{N)w^. (2.39) 

N=0 

has been used. 

The Eq.( ^.37D is useful in that all terms except the last one in the decreasing series 
do not depend on the details of any of the three string states involved (two are natural 
from the definition eq.( |2.7| )), thus knowing their levels is enough to determine them fully. 
They are hence universal. The last term Fnu is non-universal. The recursion relation 
tells that the independent functions are Xq, ■ ■ ■ , X^2 ^2, which depends on the particular 
vertex operators we inserted and so are not universal. The contribution from Fnu to F is 
generically negligible, in comparison with the universal part contribution, except the case 
A = which happens when the emitted states carry a large fraction of the total mass. 



2.2. The decay rate 

In the following, we shall calculate the decay rate ( |2.6| ). In the above discussions, we 
have fixed the levels of the incoming string states and one of the outgoing string states, 
which are Nl.r-, ^2L,2R respectively. From the mass-shell conditions, we know that once 
we fix the quanta of the incoming string states, the outgoing string states could have 
various kinds of masses, KK- momenta and windings, with respect to the energy condition 
(p75|), and conservations: 

Q_=Qi_+Q2-, Q+ = Qi+ + Q2+. (2.40) 
One important observation is that we have inequality 

+ \fWL < VnI. (2.41) 

The equality saturate when 

M Mo 

* ^ or - or- ^'-''^ 

where k is the momentum in the noncompact directions and M2, Q2- are the quanta of the 
outgoing string states with fixed level A^2L- The same inequality holds in the right-mover. 
Given a very massive initial string of high level, its state density has the asymptotic 

form 

g(N) ~ iV-^e"^ , a = 2n^^j^ . (2.43) 
The ratio between the first two terms in (|2.37] ) can be estimated to be 

g{NL -n) _ (^/NZ^-^/Nr.-2n+ml ) M) 

g{NL-2n + ml) ' ^ ' ^ 



Using the inequality (|2 . 4 1|) , it is at most of order 



Nr 



or 



exp(a(v^- V^)), v^>\^ (2-45) 

exp(a(3ViV2Z- V^)), v^<^- (2-46) 

In the extremal case N2L = Nl, one can try to calculate X directly. Thus if generically 
Nl > N2L > the first term dominates the whole summation in (|2.37| ), and the other 
terms will be neglected to get 

Fl ^ (Nl - N2L - ^ml)g{N2L). (2.47) 



Fji can be carried out in the same way. 

Note that in our approximation, Fl,r do not depend on the details of the state 
specified in eg. (p77|) and ( p.8| ) by the vertex operators, and aU states of the same level are 
emitted with the same probability. Taking advantage of this, we can get the total decay 
rate for decays into arbitrary states of given mass, winding and KK momentum by simply 
multiplying eq. (|2.6| ) by the state density ^(A^i) 



where 



and 



Remember that 



(2.48) 

AfL = NL- N2L - ^ml , Arj, = NR- N2R - ^ml, (2.49) 
. g{N^L)giN2L) . giN^R)giN2R) 

= g{N,) ' ^« = giNR) ■ ^^-''^ 

ml = ^{m^ - Ql_) ^ 2NiL. (2.51) 



As long as A^n, ^ 1 and N2L ^ 1, we can use eq.( ^.43p to write 



g, ~ (2.TH)-^(^^l^)-^e-^^^/--, (2.52) 



with the Hagedorn temperature Th = - y jyz^ and = vA^ — V-^iL — V-^2L in a 
sense coming from the kinetic energy released in the decay process. We have in the above 
restored the multiplicative constant in front of the state density. Note that for later 
convenience, in our notation we set g{M)dN = g{N)dN, a little different from usual sense 

g{M)dM = g{N)dN. 

iFrom the discussion on the inequality ( p. 41 ), we know that the dominant decay 



channel is the original incoming string break into two strings with the same 

M m M2 

and without releasing any kinetic energy k = 0. In other words, the outgoing string states 
incline to carry proportional KK-momenta and windings to their masses. Furthermore, 
in such string breaking without releasing the kinetic energy, it is least probable that the 
original string break into two strings with the same mass and quanta. 
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To see more clearly the decay process, we may take into account the dependence on 
the momentum k and other quanta. An efficient way is to expand around the maximum 
of the exponential. It turns out to be 

where e characterize the deviation of from a. Thus the dependence on k is still a 
Gaussian distribution, which could be reduced to the one in |1^] without compactification. 
The influence of the KK-momenta and windings is the change of the coefficients in the 
exponential 



expf-^./^^^^) -expf <^^J^. .2 55^ 

which tell us when a is huge, namely the original windings and KK-momentum is small, 
the distribution reduces to the one without compactification. Moreover, the contribution 
of KK-momentum and windings has also a Gaussian distribution, and the total contribu- 
tion could roughly be integrated out and give us a finite number. However, things are 
more interesting here: from the distribution of the e, it seems that when a —>■ oo, which 
happens at Q- = 0, all kinds of quanta of Q2- could give the same contribution. This is 
just an illusion, due to the subtlety on the expansion. In this case, the outgoing strings 
tends to carry very small Q2- quanta, with a distribution as exp(— "^j^^ 2!^)' where d 
characterize the deviation from the extreme. In any case, one may take into account of 
all the contributions of the KK-momenta and windings well approximately by integrating 
over the distributions above. 



With the above assumptions, eq. (|2.49|) is approximately 
NL-N2L-^ml^—-^ ^, NR-N2R--ml^'^^ ^, (2.56) 



where 3 = 

Q+ 

Now gathering all above together, we can calculate the total rate for a string with 
mass M, KK momentum n^, winding Wi to produce a string with mass m, by integrating 
over noncompact momentum k (or say it in another way, summing over the other final 
string's mass level), and summing over partitions {nu, n2i) and {wu, W2i) with constraints 
ni = flu -f n2i and Wi = wu + W2i- 
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First, we do the summations 



Tm{M,ni,Wi) = ^ ^ T[{M,ni,Wi) {m,nu,wu) + {M2,n2i,W2i)]. (2.57) 

To illustrate the picture more clearly, let us consider a special case, where the initial string 
carry very small Q-. In this case, QlQr can be divided into three parts: the kinetic energy, 
the KK modes and winding modes 

QlQr ~ Q'^g'^'^Q'', (2.58) 



where 



h"^ M 

^' = exp[-— — ], (2.59) 

21 H mM2 



expL 



2Th ^ Rf m M2 M 

i=l 



= exp[- A y + !fi _ (2.61) 

1=1 

We see from eq.( |2.59| ) that the probability to emit large kinetic energy strings is expo- 
nentially suppressed as expected. The factor ( 2.60|) from the KK modes, when we fix 



ni = nii + n2i, has a saddle point which means that the processes satisfying ^ = ^ are 
much preferred. Similarly the saddle point for the winding modes occurs at = j^. 
That is, for the final two strings, the heavier one tends to carry larger KK momentum and 
more windings. This is consistent with our general argument above. The decay rate can 
correspondingly be factorized into three parts 

T^{M,n,,w,)=T'' G'''' E (2.62) 

where we make saddle point approximation for the noncompact momentum, that is we let 
/c = 0, or M2 = M — m, when considering the KK modes and winding modes. Note that 
the basic function has the expansion 

00 

'^(^1^) exp(7rm^r + 27rmz^), (2.63) 



n= — oo 
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so the two summations in ( p.62 ) which are actuaUy summations over single argument 
Z]^i=-oo S^ii=-oo' written more concretely in terms of t? functions 



where 



A' = -?^T777?^-^L (2-66) 



Th M(M - m) 
27iiTHRf M -m'' 



^^ = 1TZ7^TT^^ (2-67) 



= ^5 ./^ . (2-68) 

_ i M 

" 27rTHRj m{M - m) ' ^^'^^^ 
, _ 2zi?2 M 
~ txTh m{M - m) ' 

In general, as we argued above, the contribution from the KK-modes and windings 
could be just a constant and the dependence on the momentum k could be encoded in 



^2 - 1/32 _ i^^M2,2 _ / , a , (3 ^ k"^ M 
72 



/92 ^ 4 ^ 'V^^^ y^^:rT^2Tj^ 2mM2 



Note that if we measure the energy of the one of the decay string with mass m, we find a 
temperature: 

Tra = ^( + r^ )-\l - -). (2.72) 

This reflects both the recoil effect and the influence of the KK-modes and windings. 

It is remarkable that when o; = 1 or /3 = 1, the decay is completely suppressed. 
This interesting limit corresponds to the BPS condition M = Q- or M = in the 
superstring case. As we show in the next section, all the discussion here could be applied 
to the superstring. Therefore, the fact that the BPS string states are stable is reflect in 
the above special limit. Actually, the suppression not only happens in the decay to the 
final massive strings, it also happens in the thermal radiation of the massless particles ||12|| . 
Moreover, from ( |2.71| ), the decay of the near BPS string states are greatly suppressed. 
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The measure of the integral in T is dN2 ~ M2(it, where we define 

f = ( _ _(_ ^ f2 731 



Expressed in an integral over t 



V ^ const. A'^X{2'kTh)~^ —m^ ^ / t — ^ — g '^dt, (2.74) 

where m^j = IHiiM^inl^ ^^d A'^^ = nf=i -^i"^- Thus including the phase space factor, we 
have a Maxwell-Boltzmann distribution for t, which could be taken as a measure of the 
total kinetic energy. Also the mean kinetic energy released per decay could be characterized 
by 

< t >= ^~^~^^ Th, (2.75) 

which is independent of the masses and satisfies the equipartition principle in D — 1 — dc 
noncompact spatial dimensions. Carrying out the integral, which is just a F-function, we 
get 

^ A'^g^-m-^^ , (2.76) 
m 

_ (D-l+dc) 

with A'j-, = const. A'^Th ^ 

Then, we take into account the fact that there are p{m)dm levels between mass m 
and m + dm, with p{m) = so the total emission rate of strings of mass m is 

(nii,n20 {'wii,W2i) 

= r V(m) = Ai^glMml-^ , (2.78) 
with A'^ = \A'j^. It can easily be checked that if we set dc = 0, the above result reduces 



where 



to the fiat space case considered in ||TJ] as expected. However, here is a puzzle. Naively, 
from semiclassical point of view, one may expect to interpret the probability for two points 
with distance m on the string to meet is inversely proportional to the volume 7tt,^-i-<^c/2^ 
where m^/^ is the mean distance between two points and D — 1 ~ dc is the noncompact 
spatial dimensions. In other words, the less the noncompact spatial dimensions, the more 
possible the two points on the string to meet, inducing the breaking of the string. But 
( p.78|) tells that this is not true: it seems that the more we compact the dimensions, the 
less possible is the string break. 
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3. Superstring decay: flat space 

In this section, we calculate the decay rate of a typical massive type II string state 
in the 10-d dimensional Minkowski space. We consider first the emission of NS-NS states, 
which we will argue later is the dominant decay channel. The basic quantity to compute 
is one piece of the semi-inclusive amplitudes squared as defined in ( |2.7|) and (pT8|), 



^i\N ^f\N 

Carry out the projection 



^ = ^E E \{^f\v{m.)f. (3.1) 



-z' 



g{Af) Jc z' Jc> z '2 



where F is the world sheet fermion number. 



z-^Tri^—^ V\k, 1)^—^ z'^V{k, l)z^], (3.2) 



We will mostly follow the procedure of and the main task is to check whether 
there is still a recursion relation when the fermions are included. 
Convert the above tree level trace to a cylinder correlator 

Tr[^^-^VHkA)^^^^z"^VikA)z'']= fMsiw) < V'HkA)^^^V'{Kv) >, 

(3.3) 

where 

with Q]^s the state density of the NS sector. Note that now the massless level already 
involves \ unit of excitation, so the factor ^= should be included. 



We define in the same way as in ( |2.18| ). On the cylinder we use vertex operators 



in the picture in the notation of |jT^ . For example, the massless state is 

Vq = {i.dX + ii4k4)e^''-^ . (3.5) 
The state generated by ip^s has the vertex operator 

2 

Vi = {i.d'^X + ii.di;k4)e'^-^ . (3.6) 

More highly excited states can be constructed with higher derivatives and more terms with 
building blocks as above. 
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These facts suggest that the correlator in ( |3.3| ) can be factorized as in ( p.24|) , imposing 
the same gauge conditions. Now the polynomial is also a function of the fermion correlators 
and their derivatives 



< V'\k,l)V'{k,v) >= i/)"^>(0,a^n...;A,a^A,...), (3.7) 

where we define 

A{v,w) =< ip{l)i;{v) > . (3.8) 

Actually there are two subtleties in arriving ( ^.Tj ) . One is that though the vertex operator 
of the massive string state might be just the linear combinations of the terms which are 
the product of chiral and anti-chiral parts, the contribution from the contraction of two 
exponentials is universal and the end result could still be written as ( ^.7|) . The other 
subtlety is the gauge choice which allow the contraction between the exponential and the 
terms without exponential vanishing. As in bosonic case ( |2.20|) , such a gauge is believed 
to exist in the critical dimension. 

What is new for the superstring case is that due to the presence of fermions, the 
correlator in ( p.3|) is a summation over spin structures. There are four spin structures 
on the cylinder. Restoring spin structures, the above correlator or Green function can be 
written as 

G{a,b)iwi - W2) =< ■ll;{wi)ljj{w2) >(a,6), (3.9) 

with a, 6 = 0, ^ labelling the periodicity 

G'(,,,)(z + 27r) =e2-«G'(„,,)(z), G'(„,fe)(z + zt) = e2-^G'(„,fe)(z). (3.10) 
The Green functions for the so-called even spin structures can be explicitly written as 

< rivJiWM (3.12) 

< ri^.m >,i.„,^ (3.13) 

We can explicitly check their periodicity by using properties of functions 

?9i(z/ + l|r) = -i?(Z/|r), (3.14) 
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it)d\{{)\it) 




zt)^?2(0|zt)' 
|zt)#i(0|it) 







M^ + l\T) = -^2il^\T), (3.15) 

^3{i^ + 1\t) = M'^\t), (3.16) 

i)^{iy + l\T) = ^^{iy\T). (3.17) 

These three spin structures are permuted under modular transformations, so to ensure 
modular invariance they should always be kept at the same time. 

For the doubly periodic case (a = 0,6 = 0), which is often called odd spin structure, 
the existence of zero modes of the Dirac operator make things subtle(see for example 
|T^ ). It is modular invariant itself and does not mix with the other three under modular 
transformation. And it is well known that this sector does not contribute to the correlation 
functions if the insertion of the vertex operator cannot eat the zero modes, while in our 
consideration of very massive strings, the corresponding vertex operators can swallow the 
zero modes without trouble. 

Note that in our convention, the arguments in the 'i?(z^|r) function read 

Inv 2ni 

v = - , T=- . 3.18 

mw inw 

To get a second contour we make the transformation v' = w~^v, which amounts to replace 
t9(z/|r) by - l|r), or G'(a,b)(^) by G'(a,6)(2 - Stt). 
Using ( |3.10| ) we get 

k{w~^v, w) = ±A(t;, w). (3.19) 

The minus sign in (|3.19| ) is irrelevant since we always have even numbers of fermions in 
our vertex operators after GSO projection. 

Taking derivatives of d functions we can see that derivatives of A get at most an extra 
irrelevant minus sign under the transformation v' = w~^v^w' = w. So it turns out that the 
fermion correlators have the same good property as dX correlators that they (effectively) 
do not change under our transformation. This is after all not a surprise since we know, for 
example from bosonization, that two fermions correspond to the derivative of a bosonic 
field, and GSO projection enables us to always have paired world-sheet fermions. 

The contour integrations on contour C'^ and around the singularity v = \ both carry 
on to the superstring case. The latter integration depends only on the asymptotic behavior 
of the correlators of e^^-^ ^ dX^ and '^'0 near f = 1, while the last two have the same order 
dependence on r 

< dX{l)dX{v) > \ ~< i(;i(;{1)iIjiIj{v) > . (3.20) 
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With proper normalization we get the same result as in the bosonic case. 

Thus we get the same recursion relation for 1^ ( |2.35| ) as in the bosonic string ; the 
only change is to replace space-time dimension 26 by 10. 

Using the expression for NS state density 

Qns{N)= j'^w-'^fNS , (3.21) 
we get the superstring semi-inclusive decay rate as a series summation 

A 1 1 

P = - m2(p - -))gNs[N -np + -m^p^ - p)] + F^u, (3.22) 

with N the mass level of the initial string; n = N — N' , N' the mass level of one of the 
final strings states; m the mass of the other of the final strings whose state is fixed from 
the start. The above expression has the same form as in the bosonic string. 

Note that we considered only the NS sector, and in the following we will argue in the 
spirit of Ijl^ that NS emission is the dominant decay channel for large, and R sector is 
suppressed by jj, with being the initial string mass level. Here for the semi- inclusive 
amplitude squared (^T|), it is more convenient to introduce only a projector for one of the 
final states 

^z'^-^' < nV{kA)^—^V{Kv)\^, > . (3.23) 
Before GSO projection, the lowest states in R sector are already massless. Their vertex 



operators are of the form [18 



k^) = e-^w"e„e*'=^, (3.24) 
where the (p term comes from superconformal ghosts, and the G term, called spin fields. 



comes from worldsheet fermion zero modes. Compare this with the massless state (|3.5| ) in 
NS sector, we see that to get to the same mass level, NS sector needs one more piece of 
excitation of the weight of dX than for R sector. The same argument works also for the 
higher mass levels. 



Relative to R sector, this extra excitation will contribute in ( p.23|) a factor of order 



< ^i\dX {l)dX {v)\^i > for each term in NS sector. Ignoring dependence on the worldsheet 
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coordinates, dXdX acts mainly as the energy-momentum tensor which can be seen by 
noting that the X scalar contribution to the energy-momentum tensor is 

T^{z) = -^dX^dX^. (3.25) 

So < (^i\dX{l)dX{y)\(^i > is of order the mass square of the initial string state which 
scales as with large. This extra factor of order makes NS sector the dominant 
decay channel. And our truncation to NS sector in the above is seen to be enough at the 
leading order. 



4. Superstring decay with compact ificat ion and KKLMMT scenario 

We can also compactify type II strings on the torus. According to the above calcula- 
tions, we will obtain the same formulas as in ( |2.37| ), ( |2.47| )-( p.50| ). Now the state density 
has the asymptotic form 

g{NL) ^ (4.1) 



The computation of the multiplicative constant in the r.h.s. is included in the appendix 
for completeness. 

The results on the toroidal compactifications can be used to estimate the massive 



string decay rates in the more realistic KKLMMT model [T^], where the type IIB strings 
live in a highly warped throat with Klebanov-Strassler geometry fl^. It was shown in 
Il20[] how to use the string reconnection probability for toroidal compact ificat ion to study 
corresponding processes in the warped geometry. We will make use of their treatment. 
A general warped geometry 

ds^ = H-^Y)7]f,^dX^dX'' + H^{Y)g,jdY'dY^ (4.2) 

will induce a potential |^ 



V{Y) = (4.3) 

27va'HHY) 

on the string worldsheet, this potential confines the string to a small region in the compact 
dimension near the position where this potential has a minimum. Fluctuations in these 



compact dimensions [2C 



■ a' 1 

< Y'Y"- >= — lull H , 1, (no sum on i) (4.4) 

2 ^ 27va'''H^0)^fV{0y ^ ^ ^ ' 
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provide an effective compactification volume for strings living in. When dfV is small in 
string units, that is, the geometry varies slowly on the string scale, we can look upon the 
strings in such a warped geometry as in a box with the effective compactification volume, 
and flat space calculations above applies. 

Near the base of Klebanov-strassler solution [|19| , the geometry is locally R'^ x with 
the warp factor 

H{Y)=H{0){1--^) (4.5) 
9s J a 



near the origin of R . H{0)~i is determined to be of order 10~ [^, and b' nearly 1 [|2l 
And J counts the number of flux. 

As in ||2^, write the minimum volume of a six-torus as Vmin = (47r^a')'^, and the 



effective volume of the six transverse dimensions as V±. The dimensionless quantity 



rmn 



V, 



3 



can be used to measure the effects of compactification. Quantum fluctuations on the R 

3 3 

contribute — t^^^^ — to ^y'" |l20[] . For S"^, quantum fluctuations give — 3^^^-*/ — , and its 



volume constrains the factor to be no less than — ^^^^ [ pO| . So flnally we get 



1^ « inf \^ (H^l! 1 (4 6) 

Vi_ (g,J)i In* (1 + g,J) ' In* (i? 4 (0)) In* (1 + g,J) 

Here we need not to worry about winding modes, and as a rough estimate we can 
also neglect the contribution from KK modes. The flat metric total decay rate ( p.48|) for 
a typical string of mass M to two strings with mass m and M2, when these two modes are 
excluded, can be simplifled 

r[M ^ m + M2] ^ Az,-d,-^||Ar26;2/,^-3-d.^ ^4 7) 

with M = Ml = A/ij, and Q = Ql = Qr, and compactiflcation volume V±_ = Yli'^ Ri- 
Remember that the initial string is averaged, and the flnal two strings both sum over given 
mass level. For the warped geometry, ([4.7] ) represents results seen by a 10- dimensional 
observer. And the quantities measured by the 4-dimensional observer, in our case the 
mass and the momentum, get red-shifted. Thus we should multiply M, m, M2 and k all 
by a red-shift factor iJ^(O). Now the mass shell conditions read like 

a'M'^H^ (0) = 4(A^ - 1). (4.8) 
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Let's look at the eq. (|2.74|) . Now we need to replace thy (0)t and tjir by if 4 (0)m/j 



to get replaced by r^if~V(0), with D = 10 and dc = Q explicitly inserted. 

Then the emission rate of a string with mass m from a typical string of mass M in 
the Klebanov-Strassler geometry can finally be estimated to be 

'Ej^^A.SS^H-'imfnn'^, (4.9) 

1 ^min 



with AkS = 16^5 v-x • 

With if- 4(0) ~ 10-4, we see that the emission rate is suppressed by a factor of order 
10-^^ in this warped geometry! With A^s no larger than 10-^, gs no larger than unity, 
the multiplicative constant in ( [4 .91) is no larger than 10-^*^. Thus in this warped geometry, 
the decay channels for a massive string to two massive strings are greatly suppressed. 

We note that the situation is different if one of the two final strings is massless. It is 



found |T3] that the decay rate is proportional to the mass of the initial string. So including 



the redshift factor, it can be written as 



Tmassless--9lMH-^{0), (4.10) 

from which we see that this decay channel is enhanced by order 10^. Thus in such warped 
geometries, massive strings may lose their energy mainly through gravitational radiation. 



5. Discussions 

We calculated the average decay rate of a typical massive string state in super- 
string theory, generalizing the earlier results obtained in the bosonic string theory. We 
did this calculation with a torus compactification, which at the first look is not realistic 
enough. However, we see that the calculation can be applied to the brane-inflation model 
of KKLMMT, we find that the decay is greatly suppressed. We believe that this shall have 
important phenomeno logical implications. 

To apply our calculation to the problems concerning cosmic strings, we need to know 
how fundamental strings are generated, for instance in the process of brane-anti-brane 
annihilation. With that knowledge, we can write down a set of differential equations on 
the distribution function of massive strings, making use of our result. We may need to take 
the evolution of the universe into account too. The approach outlined here may replace 
the usual simulation process to predict the number and the distribution of massive cosmic 
strings that remain to this day in the sky. 
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Particularly, it has been argued in [22| that gravitational radiation is not an efficient 
enough energy loss mechanism for the string network to be consistent with present obser- 
vation. Other mechanisms are needed to prevent the network from quickly dominate the 
universe. And from our calculations above, we can see that the large suppression factor is 
a generic feature for string interactions in the warped geometry with three massive strings 
involved. Thus it is important to study more carefully the energy loss mechanisms for F-D 
string networks in warped geometries. 
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Appendix A. Compute the state density for superstrings. 



The asymptotic densities of the states for bosonic and superstrings are standard text- 
book knowledge (see for example [T7] and [0). However, in our computation we need 
to know the detailed numerical coefficients for superstrings (Manes has done for bosonic 
string in H]), which are often not explicitly demonstrated, so we include this appendix to 
write them out. We follow the whole procedure of []T^ . 

Consider first open strings and NS sector only. The state degeneracy QNsi"^) is given 

by 



fNs{w) = TV 



1 + e 



i-kF 



-W 



N 



1 + w' 



(A.l) 



n=0 n=l 

with N the summation of the bosonic and fermionic number operators. Generalization of 
Hardy-Ramanujan formula gives 



n(T^)"' = ^4(0N = (-^)-t^2(0|ei^). 



(A.2) 



where the modular transformation of function 

7?4(0|r) = (-ir)-^^?2( 

has been used, with 

ilnw 

T = . 

TV 



1, 



(A.3) 



(A.4) 
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As — s> 1, the second argument of 'i?2, which now reads 

-' = -i = -T^. (A.5) 

r m w 

approaches oo. We know from the expansion 

oo 

^2(0|r') = ^ e^"("-^)'"' (A.6) 

n= — oo 

that 

i?2(0|r' ^ cx)) ^ 2e^^' = 2e^. (A.7) 
Thus ( [A .21 ) is asymptoticaUy 

r^2exp . A.8 

n=l 

^From ( |A.1|) , the state degeneracy Qns^p) can be expressed as a contour integral on 
a small circle around w = 

GNsin) = (f ^dw. (A.9) 

To compute the above integration, we make a saddle point approximation near w = 1. 
The power of w can be put on the exponential 

GNs{n) = ^ /8(-i^)42-8exp[-:^-(n + l)ln«;](i«;, (A.IO) 
2m / TT liiw 



to get the saddle point at 



where expansion can be made 



Then Qns{p) is approximately 



lnu;o = ^=, (A.ll) 



\nw = \nwQ+iu. (A. 12) 
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QNsin) ~ -^-L(^)4e-^8" / exp(-^-— (A.13) 



2tx 32 Vn y_oo TT 

Carrying out the integration over u we find 

GNsin) ~ 2-^n-^e^^^'\ (A.14) 
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Or using n ~ a'm^, write it out in terms of mass 



^7vs(m) ~2--a' ~m--e"^^"'"". (A.15) 



Here we use the convention Qis[s{m)dn = QNs{'^)dn, different from [|T7|. 

At this point, we also note that R sector has the same expression. And combine the 
left and right pieces together we arrive at the expression for closed strings 



g^\n) = [G°''{n)Y ~ 2-2n--e*"^'". (A.16) 

Taking care of the difference between the mass shell conditions of open and closed strings 
[a'rin? ~ An for closed strings) , the state degeneracy for closed string as a function of mass 
reads 

g^\m) ~ 2^a'"^m-^^e^^^"^. (A.17) 
Thus open and closed strings have the same Hagedorn temperature 

Th = (A.18) 
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